Search for: All records

ABSTRACT The dark matter (DM) distribution in dwarf galaxies provides crucial insights into both structure formation and the particle nature of DM. GraphNPE (Graph Neural Posterior Estimator), first introduced in Nguyen et al. (2023), is a novel simulation-based inference framework that combines graph neural networks and normalizing flows to infer the DM density profile from line-of-sight stellar velocities. Here, we apply GraphNPE to satellite dwarf galaxies in the FIRE-2 Latte simulation suite of Milky Way-mass haloes, testing it against both Cold and Self-Interacting DM scenarios. Our method demonstrates superior precision compared to conventional Jeans-based approaches, recovering DM density profiles to within the 95 per cent confidence level even in systems with as few as 30 tracers. Moreover, we present the first evaluation of mass modelling methods in constraining two key parameters from realistic simulations: the peak circular velocity, $$V_\mathrm{max}$$, and the peak virial mass, $$M_\mathrm{200m}^\mathrm{peak}$$. Using only line-of-sight velocities, GraphNPE can reliably recover both $$V_\mathrm{max}$$ and $$M_\mathrm{200m}^\mathrm{peak}$$ within our quoted uncertainties, including those experiencing tidal effects ($$\gtrsim 63~{{\rm per\ cent}}$$ of systems are recovered within our 68 per cent confidence intervals and $$\gtrsim 92~{{\rm per\ cent}}$$ within our 95 per cent confidence intervals). The method achieves $$10-20~{{\rm per\ cent}}$$ accuracy in $$V_\mathrm{max}$$ recovery, while $$M_\mathrm{200m}^\mathrm{peak}$$ is recovered to $$0.1-0.4 \, \mathrm{dex}$$ accuracy. This work establishes GraphNPE as a robust tool for inferring DM density profiles in dwarf galaxies, offering promising avenues for constraining DM models. The framework’s potential extends beyond this study, as it can be adapted to non-spherical and disequilibrium models, showcasing the broader utility of simulation-based inference and graph-based learning in astrophysics. 

Abstract The observed prevalence of galaxies exhibiting bursty star formation histories (SFHs) atz≳ 6 has created new challenges and opportunities for understanding their formation pathways. The degenerate effects of the efficiency and burstiness of star formation on the observed UV luminosity function are separable by galaxy clustering. However, quantifying the timescales of burstiness requires more than just the continuum UV measurements. Here we develop a flexible semi-analytic framework for modeling both the amplitude of star formation rate (SFR) variations and their temporal correlation, from which the luminosity function and clustering can be derived for SFR indicators tracing different characteristic timescales (e.g., UV continuum and Hα luminosities). Based on this framework, we study the prospect of using galaxy summary statistics to distinguish models where SFR fluctuations are prescribed by different power spectral density (PSD) forms. Using the Fisher matrix approach, we forecast the constraints on parameters in our PSD-based model that can be extracted from mock JWST observations of the UV and Hαluminosity functions and clustering bias factors atz∼ 6. If potential confusion due to e.g., dust attenuation and stellar population effects can be properly quantified, these results imply the possibility of probing the burstiness of high-zgalaxies with one-point and two-point statistics and highlight the benefits of combining long-term and short-term SFR tracers. Our flexible framework can be readily extended to characterize the SFH of high-redshift galaxies with a wider range of observational diagnostics. 

ABSTRACT The circumgalactic medium (CGM) in $$\gtrsim 10^{12}\ \mathrm{M}_{\odot }$$ haloes is dominated by a hot phase ($$T \gtrsim 10^{6}$$ K). While many models exist for the hot gas structure, there is as yet no consensus. We compare cooling flow models, in which the hot CGM flows inwards due to radiative cooling, to the CGM of $$\sim 10^{12}{\,\rm to\,}10^{13}\ \mathrm{M}_{\odot }$$ haloes in galaxy formation simulations from the Feedback in Realistic Environments (FIRE) project at $$z\sim 0$$. The simulations include realistic cosmological evolution and feedback from stars but neglect AGN feedback. At both mass scales, CGM inflows are typically dominated by the hot phase rather than by the ‘precipitation’ of cold gas. Despite being highly idealized, we find that cooling flows describe $$\sim 10^{13}\ \mathrm{M}_{\odot }$$ haloes very well, with median agreement in the density and temperature profiles of $$\sim 20{{\ \rm per\ cent}}$$ and $$\sim 10{{\ \rm per\ cent}}$$, respectively. This indicates that stellar feedback has little impact on CGM scales in those haloes. For $$\sim 10^{12}\ \mathrm{M}_{\odot }$$ haloes, the thermodynamic profiles are also accurately reproduced in the outer CGM. For some of these lower-mass haloes, cooling flows significantly overpredict the hot gas density in the inner CGM. This could be due to multidimensional angular momentum effects not well captured by our one-dimensional cooling flow models and/or to the larger cold gas fractions in these regions. Turbulence, which contributes $$\sim 10{\!-\!}40{{\ \rm per\ cent}}$$ of the total pressure, must be included to accurately reproduce the temperature profiles. Overall, cooling flows predict entropy profiles in better agreement with the FIRE simulations than other idealized models in the literature. 

Abstract Using the FIRE-2 cosmological zoom-in simulations, we investigate the temporal evolution of gas-phase metallicity radial gradients of Milky Way–mass progenitors in the redshift range of 0.4 <z< 3. We pay special attention to the occurrence of positive (i.e., inverted) metallicity gradients—where metallicity increases with galactocentric radius. This trend, contrary to the more commonly observed negative radial gradients, has been frequently seen in recent spatially resolved grism observations. The rate of occurrence of positive gradients in FIRE-2 is about ∼7% for 0.4 <z< 3 and ∼13% at higher redshifts (1.5 <z< 3), broadly consistent with observations. Moreover, we investigate the correlations among galaxy metallicity gradient, stellar mass, star formation rate (SFR), and degree of rotational support. Metallicity gradients show a strong correlation with both sSFR and the rotational-to-dispersion velocity ratio (v_c/σ), implying that starbursts and kinematic morphology of galaxies play significant roles in shaping these gradients. The FIRE-2 simulations indicate that galaxies with high sSFR ( $\log \left(\mathrm{sSFR}\right[{\mathrm{yr}}^{-1}\left]\right)\gtrsim -9.2$) and weak rotational support (v_c/σ≲ 1) are more likely—by ∼15%—to develop positive metallicity gradients. This trend is attributed to galaxy-scale gas flows driven by stellar feedback, which effectively redistribute metals within the interstellar medium. Our results support the important role of stellar feedback in governing the chemo-structural evolution and disk formation of Milky Way–mass galaxies at the cosmic noon epoch. 

ABSTRACT Galaxy formation is a complex problem that connects large-scale cosmology with small-scale astrophysics over cosmic time-scales. Hydrodynamical simulations are the most principled approach to model galaxy formation, but have large computational costs. Recently, emulation techniques based on convolutional neural networks (CNNs) have been proposed to predict baryonic properties directly from dark matter simulations. The advantage of these emulators is their ability to capture relevant correlations, but at a fraction of the computational cost compared to simulations. However, training basic CNNs over large redshift ranges is challenging, due to the increasing non-linear interplay between dark matter and baryons paired with the memory inefficiency of CNNs. This work introduces EMBER-2, an improved version of the EMBER (EMulating Baryonic EnRichment) framework, to simultaneously emulate multiple baryon channels including gas density, velocity, temperature, and H i density over a large redshift range, from $z=6$ to $z=0$. EMBER-2 incorporates a context-based styling network paired with Modulated Convolutions for fast, accurate, and memory efficient emulation capable of interpolating the entire redshift range with a single CNN. Although EMBER-2 uses fewer than 1/6 the number of trainable parameters than the previous version, the model improves in every tested summary metric including gas mass conservation and cross-correlation coefficients. The EMBER-2 framework builds the foundation to produce mock catalogues of field level data and derived summary statistics that can directly be incorporated in future analysis pipelines. We release the source code at the official website https://maurbe.github.io/ember2/. 

Abstract Integral field spectroscopy (IFS) is a powerful tool for understanding the formation of galaxies across cosmic history. We present the observing strategy and first results of MSA-3D, a novel JWST program using multi-object spectroscopy in a slit-stepping strategy to produce IFS data cubes. The program observed 43 normal star-forming galaxies at redshifts 0.5 ≲z≲ 1.5, corresponding to the epoch when spiral thin-disk galaxies of the modern Hubble sequence are thought to emerge, obtaining kiloparsec-scale maps of rest-frame optical nebular emission lines with spectral resolutionR≃ 2700. Here we describe the multiplexed slit-stepping method, which is >15 times more efficient than the NIRSpec IFS mode for our program. As an example of the data quality, we present a case study of an individual galaxy atz= 1.104 (stellar massM_*= 10^10.3M_⊙, star formation rate, SFR = 3M_⊙yr⁻¹) with prominent face-on spiral structure. We show that the galaxy exhibits a rotationally supported disk with moderate velocity dispersion ( $\sigma =36_{-4}^{+5}$km s⁻¹), a negative radial metallicity gradient (−0.020 ± 0.002 dex kpc⁻¹), a dust attenuation gradient, and an exponentially decreasing SFR density profile that closely matches the stellar continuum. These properties are characteristic of local spirals, indicating that mature galaxies are in place atz∼ 1. We also describe the customized data reduction and original cube-building software pipelines that we have developed to exploit the powerful slit-stepping technique. Our results demonstrate the ability of JWST slit-stepping to study galaxy populations at intermediate to high redshifts, with data quality similar to current surveys of thez∼ 0.1 Universe. 

ABSTRACT Fuelling star formation in large, discy galaxies requires a continuous supply of gas accreting into star-forming regions. Previously, we characterized this accretion in four Milky Way mass galaxies ($$M_{\rm halo}\sim 10^{12}{\rm M}_{\odot }$$) in the FIRE-2 cosmological zoom-in simulations. At $$z\sim 0$$, we found that gas within the inner circumgalactic medium (iCGM) approaches the disc with comparable angular momentum (AM) to the disc edge, joining in the outer half of the gaseous disc. Within the disc, gas moves inwards at velocities of $$\sim$$1–5 km s$$^{-1}$$ while fully rotationally supported. In this study, we analyse the torques that drive these flows. In all cases studied, we find that the torques in discs enable gas accreted near the disc edge to transport inwards and fuel star formation in the central few kpc. The primary sources of torque come from gravity, hydrodynamical forces, and the sub-grid $$P \mathrm{ d}V$$ work done by supernova (SN) remnants interacting with gas on $$\lesssim$$10 pc scales. These SNe remnant interactions induce negative torques within the inner disc and positive torques in the outer disc. The gas–gas gravitational, hydro, and ‘feedback’ torques transfer AM outwards to where accreting gas joins the disc, playing an important role in driving inflows and regulating disc structure. Gravitational torques from stars and dark matter provide an AM sink within the innermost regions of the disc and iCGM, respectively. Feedback torques are dominant within the disc, while gravitational and hydrodynamical torques have similar significance depending on the system/region. Torques from viscous shearing, magnetic forces, stellar winds, and radiative transfer are less significant. 

Abstract Feedback from supermassive black holes is believed to be a critical driver of the observed color bimodality of galaxies above the Milky Way mass scale. Active galactic nuclei (AGN) feedback has been modeled in many galaxy formation simulations, but most implementations have involved simplified prescriptions or a coarse-grained interstellar medium (ISM). We present the first set of Feedback In Realistic Environments (FIRE)-3 cosmological zoom-in simulations with AGN feedback evolved toz∼ 0, examining the impact of AGN feedback on a set of galaxies with halos in the mass range 10¹²–10¹³M_⊙. These simulations combine detailed stellar and ISM physics with multichannel AGN feedback including radiative feedback, mechanical outflows, and, in some simulations, cosmic rays (CRs). We find that massive (>L*) galaxies in these simulations can match local scaling relations including the stellar mass–halo mass relation and theM_BH–σrelation; in the stronger model with CRs, they also match the size–mass relation and the Faber–Jackson relation. Many of the massive galaxies in the simulations with AGN feedback have quenched star formation and elliptical morphologies, in qualitative agreement with observations. In contrast, simulations at the massive end without AGN feedback produce galaxies that are too massive and form stars too rapidly, are order-of-magnitude too compact, and have velocity dispersions well above Faber–Jackson. Despite these successes, the AGN models analyzed do not produce uniformly realistic galaxies when the feedback parameters are held constant: While the stronger model produces the most realistic massive galaxies, it tends to overquench the lower-mass galaxies. This indicates that further refinements of the AGN modeling are needed. 

Abstract Simulations and observations suggest that galaxy interactions may enhance the star formation rate (SFR) in merging galaxies. One proposed mechanism is the torque exerted on the gas and stars in the larger galaxy by the smaller galaxy. We analyze the interaction torques and star formation activity on six galaxies from the FIRE-2 simulation suite with masses comparable to the Milky Way galaxy at redshiftz= 0. We trace the halos fromz= 3.6 toz= 0, calculating the torque exerted by the nearby galaxies on the gas in the central galaxy. We calculate the correlation between the torque and the SFR across the simulations for various mass ratios. For near-equal-stellar-mass-ratio interactions in the galaxy sample, occurring betweenz= 1.2−3.6, there is a positive and statistically significant correlation between the torque from nearby galaxies on the gas of the central galaxies and the SFR. For all other samples, no statistically significant correlation is found between the torque and the SFR. Our analysis shows that some, but not all, major interactions cause starbursts in the simulated Milky Way-mass galaxies, and that most starbursts are not caused by galaxy interactions. The transition from “bursty” at high redshift (z≳ 1) to “steady” star formation state at later times is independent of the interaction history of the galaxies, and most of the interactions do not leave significant imprints on the overall trend of the star formation history of the galaxies. 

ABSTRACT We present an investigation of clustered stellar feedback in the form of superbubbles identified within 11 galaxies from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, at both cosmic noon (1 < z < 3) and in the local universe. We study the spatially resolved multiphase outflows that these supernovae drive, comparing our findings with recent theory and observations. These simulations consist of five Large Magellanic Cloud–mass galaxies and six Milky Way-mass progenitors (with a minimum baryonic particle mass of $$m_{\rm b.min} = 7100\,{\rm M}_{\odot }$$). For all galaxies, we calculate the local and galaxy-averaged mass and energy-loading factors from the identified outflows. We also characterize the multiphase morphology and properties of the identified superbubbles, including the ‘shell’ of cool ($$T\lt 10^5$$ K) gas and break out of energetic hot ($$T\gt 10^5$$ K) gas when the shell bursts. We find that these simulations, regardless of redshift, have mass-loading factors and momentum fluxes in the cool gas that largely agree with recent observations. Lastly, we also investigate how methodological choices in measuring outflows can affect loading factors for galactic winds.